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The cloning of quantum variables with continuous spectra is analyzed. A universal — or 
Gaussian — quantum cloning machine is exhibited that copies equally well the states of two conjugate 
variables such as position and momentum. It also duplicates all coherent states with a fidelity of 
2/3. More generally, the copies are shown to obey a no-cloning Heisenberg-like uncertainty relation. 
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Most of the concepts of quantum computation have 
been initially developed for discrete quantum variables, 
in particular, binary quantum variables (quantum bits). 
Recently, however, a lot of attention has been devoted 
to the study of continuous quantum variables in infor- 
mational or computation processes, as they might be in- 
trinsically easier to manipulate than their discrete coun- 
terparts. Variables with a continuous spectrum such as 
the position of a particle or the amplitude of an electro- 
magnetic field have been shown to be useful to perform 
quantum teleportation [jj] , quantum error correction , 
or, even more generally, quantum computation ||. Also, 
quantum cryptographic schemes relying on continuous 
variables have been proposed Q, while the concept of 
entanglement purification has been extended to continu- 
ous variables |5[|. In this context, a promising feature of 
quantum computation over continuous variables is that 
it can be performed in quantum optics experiments by 
manipulating squeezed states with linear optics elements 
such as beam splitters ||. 

In this Letter, the problem of copying the state of a 
system with continuous spectrum is investigated, and it 
is shown that a particular unitary transformation, called 
cloning, can be found that universally copies the position 
and momentum states. Let us first state the problem in 
physical terms. Consider, as an example of a continuous 
variable, the position x of a particle in a one-dimensional 
space, and its canonically conjugate variable p. If the 
wave function is a Dirac delta function — the particle is 
fully localized in position space, then x can be measured 
exactly, and several perfect copies of the system can be 
prepared. However, such a cloning process fails to ex- 
actly copy non-localized states, e.g., momentum states. 
Conversely, if the wave function is a plane wave with mo- 
mentum p — the particle is localized in momentum space, 
then p can be measured exactly and one can again pre- 
pare perfect copies of this plane wave. However, such a 
"plane- wave doner" is then unable to copy position states 
exactly. In short, it is impossible to copy perfectly the 
states of two conjugate variables such as position and mo- 
mentum or the quadrature amplitudes of an electromag- 
netic field. This simply illustrates the famous no-cloning 



theorem for continuous variables. 

In what follows, it is shown that a unitary cloning 
transformation can nevertheless be found that provides 
two copies of a system with a continuous spectrum, but 
at the price of a non-unity cloning fidelity. More gen- 
erally, we show that a class of cloning machines can be 
defined that yield two imperfect copies of a continuous 
variable, say x. The quality of the two copies obey a no- 
cloning uncertainty relation akin to the Heisenberg rela- 
tion, implying that the product of the x-error variance 
on the first copy times the p-error variance on the sec- 
ond one remains bounded by (fi/2) 2 — it cannot be zero. 
Within this class, a universal doner can be found that 
provides two identical copies of a continuous system with 
the same error distribution for position and momentum 
states. This doner effects Gaussian-distributed position- 
and momentum-errors on the input variable, and is the 
continuous counterpart of the universal doner for quan- 
tum bits More generally, it duplicates in a same 
manner the eigenstates of linear combinations of x and 
p, such as Gaussian wave packets or coherent states. The 
latter states are shown to be cloned with a fidelity that 
is equal to 2/3. 

In the following, we shall work in position basis, whose 
states |a;) are normalized according to (x\x') = S(x — x'). 
We assume U = 1, so that the momentum eigenstates are 
given by \p) = (27r) -1 / 2 J dxe lpx \x). Let us define the 
maximally- entangled states of two continuous variables, 
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where 1 and 2 denote the two variables, while x and p 
are two real parameters. Equation ([!]) is akin to the 
original Einstein-Podolsky- Rosen (EPR) state ||, but 
parametrized by the center-of-mass position and momen- 
tum. It is easy to check that \ip(x,p)) is maximally entan- 
gled as Tri (\ip){ip\) = Tr 2 (\ip){ip\) = 1/(2tt) for all val- 
ues of x and p, with Tri.2 denoting partial traces with re- 
spect to variable 1 and 2, respectively. The states \ip) are 
orthonormal, i.e., (tp(x' ,p')\ip(x,p)} = 8{x — x') 5{p — p'), 
and satisfy a closure relation 
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so they form an orthonormal basis of the joint Hilbcrt 
space of variables 1 and 2. Interestingly, applying some 
unitary operator on one of these two variables makes it 
possible to transform the EPR states into each other. 
Specifically, let us define a set of displacement operators 
D parametrized by x and p, 



D(x,p) = e" 
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which form a continuous Hciscnberg group. Physically, 
D(x,p) denotes a momentum shift of p followed by a po- 
sition shift of x. If D(x,p) acts on one of two entangled 
continuous variables, it is straightforward to check that 
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This will be useful to specify the errors induced by the 
continuous cloning machines considered later on. Assume 
that the input variable of a doner is initially entangled 
with another (so-called reference) variable, so that their 
joint state is \ip(0, 0)). If cloning induces, say, a position- 
shift error of x on the copy, then the joint state of the 
reference and copy variables will be \tp(x, 0)) as a result 
of Eq. (E ) . Similarly, a momentum-shift error of p will 
result in ip(0,p)). More generally, if these x and p errors 
are distributed at random according to the probability 
P(x,p), then the joint state will be the mixture 
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Let us now consider a cloning machine defined as the 
unitary transformation hi acting on three continuous vari- 
ables: the input variable (variable 2) supplemented with 
two auxiliary variables, the blank copy (variable 3) and 
an ancilla (variable 4). After applying U, variables 2 and 
3 are taken as the two outputs of the doner, while vari- 
able 4 (the ancilla) is simply traced over. Assume now 
that variable 1 denotes a reference variable that is ini- 
tially entangled with the doner input — their joint state 
is 1^(0, 0))i.2, while the auxiliary variables 3 and 4 are 
initially prepared in the entangled state 
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where f{x,p) is an (arbitrary) complex amplitude func- 
tion of x and p. The cloning transformation is defined 



-i(x4-X3) 
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where Xk (pk) is the position (momentum) operator for 
variable k. Then, the joint state of the reference, the 
two copies, and the ancilla after cloning is given by 



|*)i,|3,4 = lx ®W 2 ,3,4 |^(0, 0)>i,2 |xh,4- Using Eqs. 
and (Q), it can be written as 
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This is a very peculiar 4-variable state in that it can be 
reexpressed in a similar form by exchanging variables 2 
and 3, namely 
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with 
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g(x, P ) = —J J dx'dp' f( x ',p') (to) 

Thus, interchanging the two doner outputs amounts 
to substitute the function / with its two-dimensional 
Fourier transform g ]jTj| . This property is crucial as it 
ensures that the two copies suffer from complementary 
position and momentum errors. Indeed, using Eq. (||) 
and tracing over variables 3 and 4, we see that the joint 
state of the reference and the first output is given by 
Eq. (|J), with |/| 2 playing the role of P. Hence, the 
first copy (called copy a later on) is imperfect in the 
sense that the input variable gets a random position- 
and momentum-shift error drawn from the probability 
distribution P a (x,p) = \f(x,p)\ 2 . Similarly, tracing the 
state (^|) over variables 2 and 4 implies that the second 
copy (or copy b) is affected by a position- and momentum- 
shift error distributed as Pb(x,p) = \g(x,p)\ 2 . The com- 
plementarity between the quality of the two copies origi- 
nates from the relation between the amplitude functions 
/ and g, i. e., Eq. (|l0|), in close analogy with what was 
shown for discrete quantum doners lO] . 

Now, let us apply the cloning transformation IA on an 
input position state |xo). We simply need to project the 
reference variable onto state \xq). Applied to the initial 
joint state of the reference and the input ^(0, 0)}i,2, this 
projection operator |aro)(xo| <8> 1 yields |a;o)x|^o)2 up to a 
normalization, so the input is indeed projected onto the 
desired state. Applying this projector to the state |<&) as 
given by Eq. (§h results in the state 



dxdp f(x,p)e ipx ° \x + x) 2 \tp(x,-p)) 3A (II) 



for the remaining variables 2, 3, and 4. The state of 
copy a (or variable 2) is then obtained by tracing over 
variables 3 and 4, 
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where P a (x) — dp P a (x,p) is the position-error 
(marginal) distribution affecting copy a. Hence, the first 
copy undergoes a random position error which is dis- 
tributed as P a (x). Similarly, applying the projector on 
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the alternate expression for |<E>), Eq. (||), and tracing over 
variables 2 and 4 results in a state pb of the second copy 
that is akin to Eq. ( |l2| ) with Pb(x) = dp Pb{x,p). 
The result of cloning an input momentum state \po) 
can also be easily determined by applying a projector 
onto | — po) to the reference variable, so that the initial 
joint state of the reference and the input is projected on 
I — Po)i\Po}2- Using Eqs. (||) and (H), we obtain the 
analogous expressions for the state of copies a and 6, 



Pa(b) 



dp P a (b) (p) \Po +p)(j>o+ p\ 



(13) 



where P a (b){p) = jT^dx P a ^ b ){x,p). Consequently, 
the two copies undergo a random momentum error dis- 
tributed as P a (b){p)- The tradeoff between the quality of 
the copies can be expressed by relating the variances of 
the distributions P a (x,p) and Pb{x,p). 

Let us analyze this no-cloning complementary by ap- 
plying the Heisenberg uncertainty relation to the state 
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where [p)i(2) denote the momentum states of the first 
(second) variable. The two pairs of conjugate operators 
(xi,Pi) and (^2,^2) give rise, respectively, to the two 
no-cloning uncertainty relations 



(Ax a ) 2 (A Pb ) 2 > 1/4 
(Ax b ) 2 (A Pa ) 2 > 1/4 



(15) 



where (Ax a ) 2 and (Axb) 2 denote the variance of P a (x) 
and Pb{x), respectively. (The analogous notation holds 
for the momentum-shift distributions affecting both 
copies.) Consequently, if the cloning process induces a 
small position (momentum) error on the first copy, then 
the second copy is necessarily affected by a large momen- 
tum (position) error. 

We now focus our attention on a symmetric and uni- 
versal continuous doner that saturates the above no- 
cloning uncertainty relations. We restrict ourselves to 
solutions of the form f(x,p) = q(x) Q(—p) where Q(p) — 
(27r) -1 / 2 dx e~ l P x q(x) is the Fourier transform of 
q(x). It can be checked that this choice satisfies the sym- 
metry requirement, g(x,p) = f(x,p) Now, for the 
doner to act equally on position and momentum states, 
q(x) must be equal to its Fourier transform. Hence, the 
universal continuous doner corresponds to the choice 



f{x,p) 



1 



(16) 



so that P a ( b ){x,p) = e 



= p -(z 2 +p 2 ) 



/it is simply a bi-variate 
Gaussian of variance 1/2 on x- and p-axis. The two aux- 
iliary variables must then be prepared in the state 
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The resulting transformation effected by this universal 
doner on an input position state |sg) is given by 

\ x ) | ) |0) -» -L I r dydze-^+^Z 2 

x \x + y) \x + z) \x + y + z) (18) 

where the three variables denote the two copies and the 
ancilla, respectively. It is easy to check that Eq. ( |l8| ) im- 
plies Eq. (|12| ) and its counterpart for copy b with P a (x) = 
Pb(x) = exp(—x 2 )/y/TT, so that both copies are affected 
by a Gaussian-distributed position error of variance 1/2. 
The choice {Ax a ) 2 = {Ax b ) 2 = (A Pa ) 2 = (A Pb ) 2 = 1/2 
ensures that the doner is universal, that is, position and 
momentum states are copied with the same error vari- 
ance. The value 1/2 implies that the doner is optimal 
among the class of doners considered here in view of 
Eq. (p"5|). Furthermore, the cylindric symmetry of f(x,p) 
[i.e., it depends only on the radial coordinate (x 2 +P 2 ) 1 / 2 ] 
implies that this doner copies the eigenstates of any oper- 
ator of the form cx + dp with the same error distribution, 
as we will show. 

Let us first determine the operation of this universal 
doner on an arbitrary state |£) expressed in position ba- 
sis as J_ oo dx £(x) \x). For this, we project the refer- 
ence variable onto state i.e., the state obtained by 
changing £(x) into its complex conjugate. This projector 
I ® 1 applied on the initial state |^(0,0))i,2 yields 
the state |£*)i|£)2 up to a normalization, so the input 
is indeed projected onto Now, applying this projec- 
tor to the state |$) after cloning implies that the three 
remaining variables are projected onto the state 



dxdp f(x,p) \£(x,p)) 2 \tp{x,-p)) 3ti (19) 



where \£(x,p)) = D(x,p) |£) = f^dx' £(x')e ipx ' \x' + x) 
is the input state |£) affected by a momentum shift of p 
followed by a position shift of x. This yields 
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dx dp P a(b) {x,p) |£(x,p))(£(x,p)| (20) 



so that the two outputs arc mixtures of the \£(x,p)) 
states, with x and p distributed according to P a (h} [x,p)- 
Expressed in terms of Wigner distributions, Eq. ( p0| ) im- 
plies that W ou t(x,p) = Win(x,p)oP{x,p) witho denoting 
convolution. In particular, the universal doner simply ef- 
fects a spreading out of the input Wigner function by a 
bi-variate Gaussian of variance 1/2. 

These considerations can be easily generalized to any 
pair of canonically conjugate variables in a rotated phase 
space. First note that, using the Baker-Hausdorff for- 
mula and [x,p] = i, the displacement operator can be 
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rewritten as D(x,p) = e~ lx P/ 2 q % {p x ~ x p) _ Consider now 
any pair of observables u and v satisfying the commuta- 
tion rule [u, v] = i. Let u = cx + dp and v — —dx + cp, 
where c and d are real and satisfy c 2 + d 2 = 1. It is 
easy to check that vii — uv — px — xp, where the vari- 
ables u and v are defined just as u and v, so that D 
takes a similar form in terms of u and v (up to an irrele- 
vant phase). Therefore, as a consequence of the cylindric 
symmetry of the Gaussian |/(;r,p)| 2 , the eigenstates \u) 
of the observable u undergo a random shift of u that is 
distributed as exp(— w 2 )/y / 7r. (The position and momen- 
tum states are just two special cases of this.) We can 
also treat the cloning of coherent states (Gaussian wave 
packets) by considering the complex rotation that defines 
the annihilation and creation operators a = (x + ip)/^/2 
and a) = (x — ip)/^/2. The displacement operator can 
then be written (up to an irrelevant phase) in the usual 
form D{a) = e aa ~~ a a , where a — (x + ip)j\pl is a 
c-number that characterizes the position and momen- 
tum shift. This operator transforms the coherent state 
|ao) (i- e -i the eigenstate of a with eigenvalue ao) into 
D(a) \a Q ) = e i6 \a + a), where 9 = Im(aaJ). Thus, if 
the input of the universal doner is a coherent state |ao), 
its two outputs are a mixture of coherent states charac- 
terized by the density operator 



p = d 2 a G(a) \a + a)(a + a\ 



(21) 



where the integral is over the complex plane, and G{a) = 
2exp(— 2|a| 2 )/7r is a Gaussian distribution in a space. 
Then, using |(a|a')| 2 = exp(— \a 
calculate the cloning fidelity: 

2 



a'\ ), it is easy to 
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This fidelity does not depend on ao, so it is universal for 
all coherent states. 

Finally, we consider the cloning of quadrature squeezed 
states, defined as the eigenstates of b — (x/a + iap)/V2, 
where a is a real parameter. These states can be denoted 
as where /3 = (x/a+iap) / y/2 is a c-number. We have 

again D{p) = e^-P* 1 , so that D(J3)\/3q) = \Po+/3) up to 
a phase. In order to keep the fidelity maximum, however, 
we must use here a non-universal doner defined by 
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Both copies yielded by this doner are affected by an x- 
error of variance c 2 /2 and a p-error variance of l/(2cr 2 ), 
which implies that the density operator has the same 
form as Eq. (|l]) with G(/3) = 2 exp(-2|/3| 2 )/7r. As a 
consequence, there exists a specific cloning machine for 
each value of a that copies all squeezed states correspond- 
ing to that a with a fidelity of 2/3. In contrast, cloning 
these states using the universal doner gives a fidelity that 
decreases as squeezing increases. 



We have shown that a universal cloning machine 
for continuous quantum variables can be defined that 
transforms position (momentum) states into a Gaussian- 
distributed mixture of position (momentum) states with 
an error variance of 1/2. It is universal, as the eigen- 
states of any linear combination of x and p are copied 
with the same error distribution. In particular, it du- 
plicates all coherent states with a fidelity of 2/3. We 
conjecture that this cloning fidelity is optimal. An ex- 
perimental realization of this universal doner could be 
envisaged based on the manipulation of modes of the elec- 
tromagnetic field. The cloning transformation U would 
then couple two auxiliary modes to the input mode to be 
copied. Since IA amounts to a sequence of "continuous 
controlled-NOT gates" , it could be implemented by pair- 
wise optical QND coupling between these three modes 
As a final remark, it is worth noting that the two 
auxiliary modes must be prepared in state (|l^) , which is 
simply the product vacuum state 1 0) 3 1 0)4 processed by 
a controlled-NOT gate e~ lX3p4 . This suggests that the 
noise that inevitably arises in the cloning of the input 
mode is simply linked to the vacuum fluctuations of the 
two auxiliary modes. 

This work was supported in part by DARPA/ARO 
under grant # D A AH04-96- 1-3086 through the QUIC 
Program. N. J. C. is grateful to Samuel Braunstein, 
Jonathan Dowling, and Serge Massar for very helpful 
comments on this manuscript. 



[5] 

[6] 
[7] 

m 

[9] 
[10] 

[11] 
[12] 



S. L. Braunstein and H. J. Kimble, Phys. Rev. Lett. 80, 
869 (1998). 

S. L. Braunstein, Phys. Rev. Lett. 80, 4084 (1998). 

S. Lloyd and S. L. Braunstein, Phys. Rev. Lett. 82, 1784 

(1999). 

Y. Mu et al, Opt. Co mmun. 123, 344 (1996); T. C. 
Ralph, LAN L e-print |quant-ph/9907073| ; M. Hillery, 



LAN L e-print |quant-ph/9909006 



M. D. Reid, LANL e- 

print |quant-ph/990903C 
S. Parker, S. Bos e, and M. B. Plenio, LANL e-print 
quant-ph/990609£ . 

S. L. Braunstein, Nature (London) 394, 47 (1998). 

W. K. Wootters and W. H. Zurek, Nature (London) 299, 

802 (1982); D. Dieks, Phys. Lett. 92A, 271 (1982). 

V. Buzek and M. Hillery, Phys. Rev. A 54, 1844 (1996). 

A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. Lett. 

47, 777 (1935). 

More precisely, x — > p is an inverse Fourier transform 
while p' — > x is a direct Fourier transform. 
N. J. C erf, to appear in J Mod. Opt. (1999); also LANL 
e-print |quant-ph/9805024 . 



This solution is not unique as symmetry only requires in 
general that \g{x,p)\ 2 = \f(x,p)\ 2 . 
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